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Abstract 


This article evaluates the suitability of the ECOSSE model to estimate soil greenhouse gas (GHG) fluxes from 
short rotation coppice willow (SRC-Willow), short rotation forestry (SRF-Scots Pine) and Miscanthus after land- 
use change from conventional systems (grassland and arable). We simulate heterotrophic respiration (Ry), 
nitrous oxide (N20) and methane (CH) fluxes at four paired sites in the UK and compare them to estimates of 
R, derived from the ecosystem respiration estimated from eddy covariance (EC) and R, estimated from chamber 
(IRGA) measurements, as well as direct measurements of N?O and CH; fluxes. Significant association between 
modelled and EC-derived Rp was found under Miscanthus, with correlation coefficient (r) ranging between 0.54 
and 0.70. Association between IRGA-derived Ru, and modelled outputs was statistically significant at the Aberys- 
twyth site (r = 0.64), but not significant at the Lincolnshire site (r = 0.29). At all SRC-Willow sites, significant 
association was found between modelled and measurement-derived Ru (0.44 < r < 0.77); significant error was 
found only for the EC-derived R, at the Lincolnshire site. Significant association and no significant error were 
also found for SRF-Scots Pine and perennial grass. For the arable fields, the modelled CO, correlated well just 
with the IRGA-derived Rp at one site (r = 0.75). No bias in the model was found at any site, regardless of the 
measurement type used for the model evaluation. Across all land uses, fluxes of CH; and N-O were shown to 
represent a small proportion of the total GHG balance; these fluxes have been modelled adequately on a 
monthly time-step. This study provides confidence in using ECOSSE for predicting the impacts of future land 
use on GHG balance, at site level as well as at national level. 
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gas (GHG) emissions, has increased in recent decades 


introduction (Hastings et al., 2014). The commitment of the European 
The interest in using bioenergy crops as an alternative Union is to increase the percentage of energy from 
energy source to fossil fuels, and to reduce greenhouse renewable sources to 20% of total energy consumption 

by 2020 (EU, 2009). Under the Climate Change Act 2008 
Correspondence: M. Dondini, tel. +44 01224 273810, fax +44 01224 (Great Britain, 2008), the UK government committed to 
272703, e-mail: marta.dondini@abdn.ac.uk reduce GHG emissions by 80% in 2050 compared 
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to 1990 levels; the use of bioenergy could contribute to 
this target using dedicated ‘second generation’ (2G) 
lignocellulosic crops/plantations, including short rota- 
tion coppice (SRC), Miscanthus and short rotation for- 
estry (SRF) (Somerville et al., 2010; McKay, 2011; DECC, 
2012; Valentine et al., 2012). Consequently, a substantial 
land-use change (LUC) may occur, and it might have 
considerable environmental and economic impact (Far- 
gione et al., 2008; Searchinger et al., 2008; Gelfand et al., 
2011). 

Carbon dioxide (CO2) emissions of bioenergy had pre- 
viously been assumed to be zero (Gustavsson et al., 
1995; UK, 2008) on the assumption that emissions during 
combustion are balanced by the carbon (C) uptake dur- 
ing the growth of these bioenergy plantations, but this 
fails to take account of GHG emissions following LUC 
and subsequent crop growth. To this end, it is important 
to assess the GHG balance of bioenergy crops, particu- 
larly during the first years after conversion. 

Two approaches have been widely used to monitor 
CO, fluxes: eddy covariance (EC) and the enclosure (or 
chamber) method. Eddy covariance (McMillen, 1988; 
Aubinet et al., 2012) is a technique developed to esti- 
mate land—atmosphere exchange of gas and energy at 
ecosystem scale. The measured CO, flux, known as net 
ecosystem exchange (NEE), includes ecosystem respira- 
tion (Reco) which consists of heterotrophic (Rp) and 
autotrophic (Ra) respiration, and gross primary produc- 
tion (GPP) at ecosystem scale. As photosynthesis only 
occurs during daylight hours, the night time flux is typi- 
cally used to partition the NEE signal between GPP and 
Reco. A flux-partitioning algorithm that defines a short- 
term temperature sensitivity of Reco is applied to extrap- 
olate CO, fluxes from night to day (Reichstein et al., 
2005). In a plant removal experiment (Hardie et al., 
2009), the total Rp from the whole soil profile was found 
to be approximately between 46 and 59% of the total 
Reco. Abdalla et al. (2014) used these values to simulate 
Ry, from selected European peatland sites using a soil 
process-based model, ECOSSE. 

Enclosure methods have been developed to measure 
CO, efflux from soil; these methods involve covering an 
area of soil surface with a chamber and the soil CO, 
efflux can be determined using two main modes: 
dynamic (closed or open) and closed static. In the for- 
mer mode, a steady stream of air is pumped directly in 
to the chamber (Christensen, 1983; Skiba et al., 1992). 
The latter mode simply involves closing the chamber 
for approximately 20-60 min and taking gas samples at 
intervals for analysis (Hutchinson & Mosier, 1981), or 
circulating the chamber air through a nondestructive 
infrared gas analyser (IRGA) for approximately 2 min 
(Norman et al., 1992; Smith & Mullins, 2000). Several 
studies have used the closed chamber method com- 


bined with root-exclusion methods, tree grilling or 
stable isotopes to understand the relative contribution 
of Ra and Ra to total soil respiration (Riot) under differ- 
ent land uses. 

Byrne & Kiely (2006) demonstrated that R, under 
grassland soil in Ireland accounted for approximately 
50% of Rio during the summer months and 38% during 
the rest of the year. Pacaldo et al. (2013) reported a con- 
tribution of Ra of about 18-33% of Riot under SRC-Wil- 
low at three different development stages in the USA. 
In a study on commercial farms located across the UK, 
Koerber et al. (2010) reported a contribution of Ra on 
Riot for wheat of approximately 32% from January to 
May, 79% from June to September and 67% from Octo- 
ber to December. A meta-analysis of soil respiration 
partitioning studies reported values for the ratio Ry/Rtot 
for forest soils as ranging from 0.03 to 1.0 (Subke et al., 
2006). Overall, the ratio was higher for boreal coniferous 
forests than temperate sites. In temperate, mixed decid- 
uous forests ranges for Ru / Riot of 0.3-0.6 were reported 
(Gaudinski et al., 2000; Borken et al., 2006; Millard et al., 
2010; Heinemeyer et al., 2012). Several studies have also 
shown that bioenergy plantations have low nitrous 
oxide (N2O) emissions compared to agricultural crops 
because of their lower nutrient requirements, thus 
reducing the fertilizer requirements, and more efficient 
nutrient uptake, thus increasing competition with 
microbial organisms of NO production (Flessa et al., 
1998; Hellebrand et al., 2010; Drewer et al., 2012). 

Methane (CH4) is another important GHG that may 
be a substantial component of the GHG balance from 
several terrestrial ecosystems (van den Pol-van Dasse- 
laar et al., 1999). In agricultural systems, soil is typically 
a small net source or sink for CH, (Boeckx & Van 
Cleemput, 2001). Bioenergy crops usually present either 
a small CH; sink (Hellebrand et al., 2003; Kern et al., 
2012) or a small CH, source (Gelfand et al., 2011). The 
magnitude of the CH; flux is typically much smaller 
than CO, and N20, in both agricultural soils (Boeckx & 
Van Cleemput, 2001) and bioenergy crops (Hellebrand 
et al., 2003). However, very few studies (Hellebrand 
et al., 2003; Gelfand et al., 2011; Kern et al., 2012) have 
reported on the contribution of CH, emission from 
bioenergy systems, increasing uncertainty in the direc- 
tion of this small flux (Zona et al., 2013). 

Several factors control the GHG emissions of both 
bioenergy and conventional crops, such as site manage- 
ment, for example fertilization (Crutzen et al., 2008; 
Hellebrand et al., 2008, 2010), previous land use (Smith 
& Conen, 2004) and climatic conditions (Flessa ef al., 
1998; Hellebrand et al., 2003). Despite the high variabil- 
ity of the GHG fluxes, to our knowledge, only one study 
in the UK (Drewer et al., 2012) has reported on all three 
GHG fluxes (COs, NO and CHA from soils under 
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bioenergy crops (Miscanthus and SRC-Willow) and, in 
particular, after transition from former conventional 
systems. To fill this gap, soil models are a useful tool to 
predict GHG fluxes when site measurements are not 
available, especially when studying the effects of the 
change in land use over time and under different cli- 
matic conditions over large areas. 

However, soil models need to be extensively tested 
under a range of climates and soils before being applied 
under conditions different from those used to parame- 
terize and calibrate the model itself. In fact, model eval- 
uation involves running a model using input values 
that have not been used during the calibration process, 
demonstrating that it is capable of making accurate sim- 
ulations under a wide range of conditions (Moriasi 
et al., 2007). A model can only be properly evaluated 
against independent data and a useful model should be 
able to simulate those data with some degree of accu- 
racy (Smith & Smith, 2007). 

Although several soil models have been developed 
for conventional agricultural and forest systems, most 
of them have not been fully parameterized and effec- 
tively tested for application on 2G bioenergy crops, such 
as Miscanthus, SRF and SRC (Dimitriou et al., 2012; 
Borzecka-Walker et al., 2013; Robertson et al., 2015). 
Here, we focus on the applicability of the process-based 
model ECOSSE to predict soil CO, (heterotrophic respi- 
ration), N-O and CH, after transition from conventional 
to bioenergy crops. 

The ECOSSE model was developed mainly to simulate 
the C and nitrogen (N) cycles using minimal input data 
on both mineral and organic soils (Smith et al., 2010a,b). 
The ECOSSE model has been previously evaluated across 
the UK to simulate the effect on soil C of LUC to SRF 
(Dondini et al., 2015a), Miscanthus and SRC-Willow (Don- 
dini et al., 2015b), to simulate soil NO emissions in crop- 
land sites in Europe (Smith et al., 2010b; Bell et al., 2012) 
and CO, emissions from peatlands (Abdalla et al., 2014). 

This article evaluates the suitability of ECOSSE for 
estimating soil GHG fluxes from SRC-Willow, SRF-Scots 
Pine and Miscanthus soils in the UK after LUC from 
conventional systems (grassland and arable). Based on 
previously published recommendations, a combination 
of graphical techniques and error statistics has been 
used for model evaluation (Moriasi et al., 2007). Model 
testing is often limited by the lack of field data to which 
the simulations can be compared (Desjardins et al., 
2010). In this study, the model is evaluated against 
2 years of observations at four locations in the UK, com- 
prising one transition to SRF-Scots Pine, three transi- 
tions to SRC-Willow and two transitions to Miscanthus. 
Modelled GHG fluxes from conventional systems have 
also been evaluated against field measurements (three 
grassland and two arable fields). 


Materials and methods 


ECOSSE model 


The ECOSSE model includes five pools of soil organic matter, 
each decomposing with a specific rate constant except for the 
inert organic matter (OM) which is not affected by decomposi- 
tion. Decomposition is sensitive to temperature, soil moisture 
and vegetation cover; soil texture (sand, silt and clay), pH and 
bulk density of the soil along with monthly climate and land- 
use data are the inputs to the model (Coleman & Jenkinson, 
1996; Smith et al., 1997). The ECOSSE model is able to simulate 
C and N cycle for six land-use categories of vegetation: arable, 
grassland, forestry, seminatural, Miscanthus and short rotation 
coppice willow (SRC-Willow). 

The vegetation input to the soil (SI) is estimated by a subrou- 
tine in the ECOSSE model which uses a modification of the 
Miami model (Lieth, 1972), a simple model that links the cli- 
matic net primary production of biomass (NPP) to annual mean 
temperature and total precipitation (Grieser et al., 2006). For a 
full description of the ECOSSE model and the plant input, esti- 
mates refer to Smith et al. (2010a) and Dondini et al. (2015b). 

The minimum ECOSSE input requirements for site-specific 
simulations are as follows: 


Climate /atmospheric data: 
e 30-year average monthly rainfall, potential evapotranspira- 
tion (PET) and temperature, 


e Monthly rainfall, temperature and PET. 
Soil data: 


e Initial soil C content (kg ha”), 

e Soil sand, silt and clay content (%), 
e Soil bulk density (g om "1 

e Soil pH and 

e Soil depth (cm). 


Land-use data: 
e Land use for each simulation year. 


The initialization of the model is based on the assumption that 
the soil column is at steady state under the initial land use at the 
start of the simulation. Previous work has used soil organic car- 
bon (SOC) measured at steady state to determine the plant inputs 
that would be required to achieve an equivalent simulated value 
(e.g. Smith et al., 2010a). This approach iteratively adjusts plant 
inputs until measured and simulated values of SOC converge. In 
the absence of additional measurements, estimated plant inputs 
were calculated from a feature built in the ECOSSE model which 
combine the NPP model Miami (Lieth, 1972, 1973), land-manage- 
ment practices of the initial land use and measured above-ground 
biomass (details are given in Dondini et al., 2015b). 


Data 


In 2011-2013, four sites were sampled in Britain using a paired 
site comparison approach (Keith et al., 2015; Rowe et al., 2015). 
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The sites and the relative measurements contribute to the 
ELUM (Ecosystem Land Use Modelling & Soil Carbon GHG 
Flux Trial) project (Harris et al., 2014). Each site consisted of 
one reference field (arable or grassland, depending on the pre- 
vious land use of the bioenergy fields) and one or more adja- 
cent bioenergy fields (Miscanthus, SRC-Willow, SRF-Scots Pine), 
for a total of six transitions to bioenergy at four site across UK 
(Table 1). A full description of the sites can be found in Drewer 
et al. (2012, 2015); J. McCalmont, N. McNamara, I. Donnison 
and J. Clifton-Brown (in preparation); and Z. M. Harris, G. 
Alberti, J. R. Jenkins, E. Clark, R. Marshall, R. Rowe, N. McNa- 
mara and G. Taylor (in preparation). 

At each bioenergy and reference field, the NEE data were 
obtained from continuous EC measurements (McMillen, 1988; 
Aubinet et al., 2012) using open path IRGAs (LI-7500) and sonic 
anemometers. All details regarding the EC data corrections, 
quality control, footprint and gap filling procedures can be 
found in Aubinet et al. (2003). The night time fluxes were used 
to partition the NEE flux measurements into GPP and Reco 
(Reichstein et al., 2005). 

Soil GHG fluxes were measured on a monthly basis at eight 
points randomly distributed within each field. Soil CO, fluxes 
were measured using an IRGA connected to an SRC-1 soil 
respiration chamber (PP Systems, Amesbury, MA, USA). Mea- 
surements of soil CH, and NO fluxes were made using a static 
chamber method (approx. 30 1) with the addition of a vent to 
compensate for pressure changes within the chamber during 
times of sampling. Gas samples were analysed by gas chromato- 
graph. All details regarding the chamber data can be found in 
Drewer et al. (2012), Yamulki et al. (2013) and Case et al. (2014). 

Measurements of soil C, soil bulk density and soil pH to 1 m 
soil depth, as well as information on the land-use history, were 
collected for each field (Keith et al., 2015; Rowe et al., 2015). Soil 
texture was measured for each site up to a depth of 30 cm; val- 
ues to 1 m soil depth were extracted from the soil database 
(1 km resolution) described in Bradley et al. (2005), which is a 
collated soils data set for England and Wales, Scotland and 
Northern Ireland. Air temperature and precipitation data at 


each location were extracted from the E-OBS gridded data set 
from the EU-FP6 project ENSEMBLES, provided by the 
ECA&D project (Haylock et al., 2008). This data set is known as 
E-OBS and is publicly available (http://eca.knmi.nl/). For each 
location, monthly air temperature and precipitation for the 
30 years before measurements started were used to calculate a 
long-term average (Table 2). At each site, air temperature and 
precipitation were collected during the entire study period and 
monthly values were used as input to the model. Monthly PET 
was estimated using the Thornthwaite method (Thornthwaite, 
1948), which has been used in other modelling studies when 
direct observational data have not been available (e.g. Smith 
et al., 2005; Dondini et al., 2015a). 


Model evaluation and statistical analysis 


Monthly simulations of soil CO2, NO and CH; fluxes were 
evaluated against monthly chamber measurements. In addition, 
the soil CO, predicted by the ECOSSE model was compared to 
estimates of Rp derived from the NEE measured by the EC. 

At each site, the ECOSSE model has been run for the refer- 
ence field (i.e. no land-use transition) and the bioenergy crop 
field (i.e. following transition from the reference land cover). 
The reference fields have been run for the conventional crop 
(arable, grassland) with no LUC, and the length of the simula- 
tions has been defined by the age of the plantation. At the 
bioenergy sites, the model has been run for the reference fields 
(conventional crop) with LUC to bioenergy crop; the length of 
the simulations was based on the time after transition to bioen- 
ergy crop. Measured soil characteristics and meteorological 
data have been used as inputs to drive the model (see above 
for input details), and the results of the simulations were 
compared to the GHG fluxes measured at the sites. 

We expected a monthly underestimate of the soil CO, flux 
simulations because the ECOSSE model simulates Ru (from 
living micro-organisms + decomposition of old C sources, i.e. 
saprotrophic), while the CO, fluxes measured at the sites repre- 
sent the total CO, efflux from the soil profile (Ra + Rp, chamber 


Table 1 Details of soil C, soil bulk density and soil pH to 1 m soil depth, as well as information on the land-use history at the study 
fields. Soil texture to 1 m soil depth was extracted from the soil database (1 km resolution) described in Bradley et al. (2005) 


Latitude, Establishment Bulk density 

Site Land use longitude year Carbon (%) Nitrogen (%) (g om " 
West Sussex Short rotation 50.9,—0.4 2008 0.63 0.17 1.50 

coppice (SRC)-Willow 

Grassland 50.9, —0.4 2000 0.53 0.17 1.55 
East Grange Short rotation forestry 56.0, —3.6 2009 0.95 0.18 1.47 

(SRF)-Scots Pine 

Grassland 56.0, —3.6 2009 1.30 0.17 1.49 

SRC-Willow 56.0, —3.6 2009 1.57 0.17 1.38 

Arable 56.0, —3.6 Pre-1990 1337 0.18 1.57 
Lincolnshire SRC-Willow 53.1, —0.3 2006 1.26 0.11 1.41 

Miscanthus 53.1, —0.4 2006 1.30 0.13 1.53 

Arable 53.1, —0.5 Pre-1990 1.47 0.13 1:37 
Aberystwyth Miscanthus 52.4, —4.0 2012 0.98 0.25 1.21 

Grassland 52.4, —4.0 Pre-2007 1.16 0.26 1.45 
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Table 2 Long-term (30 years) monthly rainfall, temperature, potential evapotranspiration (PET). Monthly rainfall and temperature 
were extracted from the E-OBS data set (Haylock et al., 2008; http:/ /eca.knmi.nl/). Monthly PET was estimated using the Thornth- 


waite method (Thornthwaite, 1948) 


Aberystwyth East Grange Lincoln West Sussex 

Rain Temperature PET Rain Temperature PET Rain Temperature PET Rain Temperature PET 
Month (mm) (C) (mm) (mm) (°C) (mm) (mm) (°C) (mm) (mm) (C) (mm) 
January 152 4 15 103 3 11 48 4 13 80 5 16 
February 112 4 17 72 3 15 37 4 17 54 5 18 
March 124 5 29 74 5 27 4l 6 30 55 7 30 
April 86 7 45 53 7 47 43 9 48 46 9 48 
May 82 10 69 61 10 72 45 12 73 47 12 73 
June 93 13 89 60 13 96 56 14 97 48 15 95 
July 105 15 101 67 14 105 49 17 112 49 17 110 
August 114 14 93 77 14 96 55 17 103 52 17 103 
September 121 13 71 84 12 70 49 14 76 60 15 79 
October 174 10 46 100 9 43 55 11 46 99 12 51 
November 171 7 27 94 5 22 53 7 25 88 8 29 
December 168 4 17 91 3 12 51 4 14 86 6 18 


measurements) or NEE (EC measurements). To compare the 
modelled and measured Ry, we estimated the Rp as a propor- 
tion of the measured CO) flux, depending on the measurement 
type (except EC data), vegetation type and growing season. 

The EC measurements of NEE were used to derive Reco; to 
our knowledge, only the study by Abdalla et al. (2014) has 
reported estimates of Ry from Reco. Abdalla et al. (2014) applied 
the approach proposed by Hardie et al. (2009) for peaty soils 
and reported a contribution of Rp to Reco of 46-59%. 

To represent the variations in R, throughout the year, 
Abdalla et al. (2014) assumed that R, was at the lowest value of 
the range (46% Reco) during the summer (June-August), the 
highest value (59% Reco) during the winter (December—Febru- 
ary) and at the mean value (52.5% Reco) during the rest of the 
year (March May and September November), In this study, 
we used the same approach of Abdalla et al. (2014) to derive Ry 
from EC measurements from all land-use systems. 

Chamber measurements represent the total CO, flux from 
the soil as the sum of Ra and Rp, with the exception of grass- 
land where exclusion of full leaves from the chamber is diffi- 
cult, and therefore, above-ground plant respiration is also 
included in the measurements. We conducted a literature 
review to determine the partitioning of Riot measured by the 
chambers under different vegetation types. Additional experi- 
ments within the ELUM project were also undertaken to 
directly quantify Rp and R, at selected network sites (data not 
shown); where available, we used the R, site data to estimate 
Ry from Riot measured by the chambers (Lincolnshire — Mis- 
canthus, West Sussex — SRC-Willow, Aberystwyth — Miscant- 
hus). An overview of the data source and the monthly 
proportion of Rn for each vegetation type and at each site are 
shown in Table 3. 

A quantitative statistical analysis was undertaken to deter- 
mine the coincidence and association between measured and 
modelled values, following methods described in Smith et al. 
(1997) and Smith & Smith (2007). The statistical significance of 
the difference between model outputs and experimental 


observations can be quantified if the standard error of the mea- 
sured values is known (Hastings et al., 2010). The standard 
errors (data not shown) and 95% confidence intervals around 
the mean measurements were calculated for all field sites. 

The degree of association between modelled and measured 
values was determined using the correlation coefficient (r). 
Values for r range from —1 to +1. Values close to —1 indicate a 
negative correlation between simulations and measurements, 
values of 0 indicate no correlation and values close to +1 indi- 
cate a positive correlation (Smith & Smith, 2007). The signifi- 
cance of the between simulations and 
measurements was assigned using a Student’s t-test as outlined 
in Smith & Smith (2007). 

Analysis of coincidence was undertaken to establish how 
different the measured and modelled values were. The degree 
of coincidence between the modelled and measured values was 
determined using the lack of fit statistic (LOFIT), and its signifi- 
cance was assessed using an F-test (Whitmore, 1991) indicating 
whether the difference in the paired values of the two data sets 
is significant. The EC measurements were not replicated, so the 
coincidence between measured and modelled values was deter- 
mined using the mean difference (M), calculated as the sum of 
the differences between measured and modelled values and 
divided by the total number of measurements (Smith et al., 
1997). The variation across the different measurements was 
then used to calculate the value of Student’s t-test and com- 
pared to the ¢ distributions (two-tailed test) to obtain the proba- 
bility that the mean difference is statistically significant. All 


association 


statistical results were considered to be statistically significant 
at P < 0.05. 


Results 


The ECOSSE model was evaluated by comparing the 
outputs to the EC-derived and IRGA-derived R, fluxes 
from eleven fields over four sites, representing the 
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Table 3 Contribution of heterotrophic respiration (Rp) on total respiration (Riot) at the study sites 


Arable SRC-Willow Grassland SRF-Scots Pine 
Koerber Pacaldo et al. Byrne & 
et al. (2010) (2013) Miscanthus Kiely (2006) Millard et al. (2010) 
Lincolnshire January 32% Riot 75% Riot 41% Riot 
February 32% Riot 75% Riot 41% Riot 
March 32% Riot 75% Riot 85% Riot 
April 32% Riot 75% Riot 85% Riot 
May 32% Riot 75% Riot 85% Riot 
June 79% Riot 75% Riot 85% Riot 
July 79% Riot 75% Riot 44% Riot 
August 79% Riot 75% Riot 44% Riot 
September 79% Riot 75% Riot 44% Riot 
October 67% Riot 75% Riot 44% Riot 
November 67% Riot 75% Riot 41% Riot 
December 67% Riot 75% Riot 41% Rot 
West Sussex January 82% Riot 60% Riot 
February 82% Riot 60% Riot 
March 82% Riot 60% Riot! 
April 82% Riot 60% Riot 
May 82% Riot 60% Riot 
June 82% Riot 40% Rf 
July 82% Riot 40% Riot! 
August 82% Riot 40% Riot! 
September 82% Riot 60% Riot 
October 82% Riot 60% Riot 
November 82% Riot 60% Riot 
December 82% Riot 60% Riot | 
Aberystwyth January 62% Riot 60% Riot 
February 62% Riot 60% Riot 
March 36% Riot 60% Riot! 
April 36% Riot 60% Riot 
May 36% Riot 60% Riot 
June 36% Riot 40% Riot! 
July 36% Riot 40% Riot! 
August 36% Riot 40% Riot! 
September 36% Riot 60% Riot 
October 36% Riot 60% Riot 
November 62% Riot 60% Riot 
December 62% Ru 60% Riot 
East Grange January 32% Riot 25% Riot 60% Riot 61% Riot 
February 32% Riot 25% Riot 60% Riot 61% Riot 
March 32% Riot 25% Riot 60% Riot 61% Riot 
April 32% Riot 25% Riot 60% Riot 61% Riot 
May 32% Riot 25% Riot 60% Riot! 61% Riot 
June 79% Riot 25% Riot 40% Riot! 61% Riot 
July 79% Riot 25% Riot 40% Riot! 61% Riot 
August 79% Riot 25% Riot 40% Riot! 61% Riot 
September 79% Riot 25% Riot 60% Riot 61% Riot 
October 67% Riot 25% Riot 60% Riot 61% Riot 
November 67% Riot 25% Riot 60% Riot! 61% Riot 
December 67% Riot 25% Riot 60% Riot 61% Riot 


*Values derived from direct measurements on root-exclusion plots. 
+Where Riot is 60% of measured CO, to account for plant respiration. 
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following land-use systems: grassland (permanent), ara- 
ble (barley), Miscanthus, SRC-Willow and SRF-Scots 
Pine. 

Soil CO, fluxes under Miscanthus were measured at 
two sites, Lincolnshire and Aberystwyth. At both sites, 
the modelled R, followed the same seasonal pattern of 
measured data (Fig. 1). At the Lincolnshire site, a statis- 
tically significant association between modelled and EC- 
derived Rp (r = 0.54) was found, but a small significant 
bias in the model simulations when tested against the 
EC-derived Ru was also found (Table 4). On the other 
hand, the IRGA-derived Rp did not correlate well with 
the modelled outputs (r = 0.29), but no bias was found 
in the model simulations (Table 4). 

At the Aberystwyth site, significant association 
between modelled and measurement-derived Rp, was 
found, regardless the type of measurement used. A 
slightly higher correlation coefficient was calculated 
correlating the modelled Rj, with the EC-derived Rp 
(r = 0.70) compared to the one arising from the corre- 
lation with the IRGA-derived Rp (r = 0.64). No signifi- 
cant error between simulated and IRGA-derived Rp 


was found for this site, but a bias in the model was 
found when it was tested against the EC-derived Rp 
(Table 4). 

The model performance to simulate soil CO, fluxes 
under SRC-Willow was tested against measurements 
taken at three sites: Lincolnshire, West Sussex and East 
Grange (Fig. 2). At all sites, a good agreement was 
found between simulations and measurement-derived 
Ry with r values ranging from 0.44 to 0.77. Also, no sig- 
nificant error between simulated and measurement- 
derived Rp was found, with the exception of the EC- 
derived R, at the Lincolnshire site (Table 4). 

Model performance to simulate soil CO; fluxes under 
SRF-Scots Pine has been evaluated against data collected 
at the East Grange site (Fig. 3). The modelled outputs fol- 
lowed the same pattern of the measured values, and the 
statistical analysis showed good correlation with both 
IRGA- and EC-derived Ru. Moreover, we found no statis- 
tically significant error between modelled and measured 
values as well as no bias in the model (Table 4). 

Model simulations of soil Rj, have also been evaluated 
for conventional crops (arable and grassland). Overall, 
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Fig. 1 Eddy covariance derived (dotted line with diamond markers), IRGA derived (filled triangle) and modelled (solid line with circle 
markers) monthly heterotrophic CO» (Rp) under Miscanthus plantations during the measurement period. 
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Fig. 2 Eddy covariance derived (dotted line with diamond markers), IRGA derived (filled triangle) and modelled (solid line with cir- 
cle markers) monthly heterotrophic CO; (Ry) under SRC-Willow plantations during the measurement period. 


the simulated CO; follows the same pattern as the mea- 
sured values at all sites (Figs 4 and 5). The statistics 
highlighted a significant correlation (ranging between 
0.48 and 0.87 across all sites and measurements types) 
and no significant error between modelled and mea- 
sured values as well as no model bias under perennial 
grass (Table 4). For the arable fields, the modelled CO, 
was significantly correlated to the measured value just 
for the IRGA-derived R, at the Lincolnshire site 


(r = 0.75); however, no bias in the model was found at 
any site, regardless of the measurement types used for 
the model evaluation (Table 4). 

Monthly fluxes of CH, and N2O were shown to be 
highly variable, both spatially and temporally, across all 
land uses, so we present an example of the correlation 
between modelled and measured soil N-O and CH, 
fluxes for each land use. Both N-O and CH; are very 
small fluxes and the model outputs were within the 
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Fig. 3 Eddy covariance derived (dotted line with diamond markers), IRGA derived (filled triangle) and modelled (solid line with 
circle markers) monthly heterotrophic CO (Rn) under short rotation forestry-Scots Pine plantation during the measurement period. 
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Fig. 4 Eddy covariance derived (dotted line with diamond markers), IRGA derived (filled triangle) and modelled (solid line with 
circle markers) monthly heterotrophic CO» (Rn) under arable plantations during the measurement period. 
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Fig. 5 Eddy covariance derived (dotted line with diamond markers), IRGA derived (filled triangle) and modelled (solid line with circle 
markers) monthly heterotrophic CO, (Rp) under grassland plantation during the measurement period. 


errors of the measurements, for both GHGs and at all 
sites (data not shown). low correlation 
between measured and modelled values has been found 
for the majority of the sites, ranging from —0.02 to 0.61 
for N20 and from —0.29 to 0.53 for CH4. The high vari- 
ability of the measured N-O and CH; fluxes led to a 
statistically significant error between simulated and 
measured values at most of the study sites (Tables 5 
and 6). 


However, 


Discussion 


Soil CO, emissions under Miscanthus have been quanti- 
fied at two sites (Lincolnshire and Aberystwyth) using 
two different sampling methods (EC and IRGA meth- 
ods). At both sites, we found a high correlation 
between measured and modelled Ry, ranging from 0.54 
to 0.60, except for the IRGA values at Lincolnshire site 
(r = 0.29, Table 4). The lack of association at this site 
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% was mainly due to differences between modelled and 
zis g S qs IRGA-derived Ry in the year 2013 (Fig. 1b). In April 
rs S 2013, the soil was harrowed and disked to break up 

S the rhizomes for improved yield, so the system was 

BS out of balance; the farmer also applied waste wood 

D E = a Se products, which led to high CO, emissions, undetected 
2 = A E & = ZZ by the model (May August 2013 in Fig. 1b) as this was 
< ajl È not included in the management file. In the ECOSSE 
model, the patterns of C and N debris return during 

7 êb d 8s 8 S3 = the growing season follow a standard exponential rela- 

3 & [2 3 = a > tionship, as originally derived by Bradbury etal. 
KR (1993). Any alteration, such as harrowing or waste 

D g application, cannot be easily entered by the user. The 

E E scope of the present study is to evaluate the model 

E SNE Se E using independent data which has not been used to 

a/2 2 g S Sa S develop the model. Therefore, we deliberately chose 

KR À T not to apply any modifications to the model to fit the 

X Se measured data. However, the model was able to simu- 

A 22/8 a 8 aa £ late independent data derived from two different 
fea|= =H À °= N S sources with a good degree of accuracy. 
E Soil CO, emissions under SRC-Willow and SRF-Scots 

SL Ki Pine plantations have been quantified using the same 

$ E g sampling methods. At all sites, the modelled Rp signifi- 
Sp 2 yla 8 8 eg a cantly correlated with all types of measurements, show- 
PE GIE = ia FTN g ing no significant error between measured and 
E modelled values (Fig. 2). 

e S o A 2 SE? E The model has also been tested against CO, fluxes 
z A Sg, S= 3 measured under conventional crops. At all three grass- 
£ land sites (West Sussex, Aberystwyth and East Grange), 
2 v, P the measured CO, fluxes correlate significantly with the 
T z Elg 8 9 gq Y modelled values and the statistical analysis showed no 
F GEI, S A Glatz? 5 error between measured and modelled values, and no 
2 S bias in the model (Fig. 5). This is a striking result which 
= z £ = underlines the good quality of the data provided for the 
N a z S model evaluation, as well as the good model perfor- 

ž J gje g 8 8s g mance to simulate soil CO, fluxes. 
S S Die sg SN aS 5 Under grassland, R, derived from the IRGA measure- 
Z m ments does not always show a high correlation with the 
E g S modelled values, particularly during the summer 
4 z e months (Fig. 5). This lack of correlation is mainly due to 
S els gf 2 to a the difficulties in the separation of soil respiration from 
“ 5 F Ke N m Te 2 grassland, due to the possible inclusion of vegetation 
g E within the chamber. When deriving Rn from grassland, 
S S $ we estimated that 60% of the measured CO, can be 
£ E F A attributed to plant (leaf) respiration, as reported by 

= 5 : ! 
Z 5 S a pn B D g E Byrne & Kiely (2006), but this crude estimate does not 
e Š EIERE SE : always reflect the field conditions. For an accurate quan- 
Z £ tification of the proportion of the CO, derived from the 
m e e S plant occluded in the chambers, field experiments 
6 RH 5 = x KS S would be needed to explicitly quantify plant respiration 
S S a $ SS 5 KR Kl g and biomass. 

e S E i E 3 S e S 5 7 E The analysis of the soil Ra fluxes from the arable 
ra v lsa 4722 5 28! S fields reveals reasonable model performance at the Lin- 
d S F 4 Sle g Saw 2 S < colnshire site, while at the East Grange site, correlation 
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g between modelled and measured IRGA values was poor 
2 S 23 A & SSe (Table 4). This discrepancy between modelled and 
eS et e SE measurement-derived Rp appears to be due to the nat- 
ure of the source data; in fact, the IRGA-derived Rp is 
E estimated from a single data point which is taken to 
v £ = 8 = a represent monthly CO, fluxes. Therefore, the monthly 
E S = x Ee Sea CO, flux might not be properly represented if high flux 
< Aj | variation occurred within the month. Another explana- 
tion could also be the discontinuity of the IRGA mea- 
S êp dë surements taken at the East Grange site (Fig. 4b). The 
A & S 3 z S © o E latter hypothesis is supported by the Ru results of the 
R arable field at the Lincolnshire site. In fact, the IRGA 
BS g measurements at the Lincolnshire site have been taken 
È 5 over a 2-year period, and the statistical analysis shows a 
g E good correlation against the model output (r = 0.75; 
g z z D ER S Table 4). Therefore, we conclude that the low correla- 
SE ci c sua F tion at the East Grange arable field is mainly due to the 
y variability and quantity of the measurements, and that 
gu S £ the model accurately describes the CO2 emissions from 
TSSA = 8 BB. E ble crop. 
GEAlS a a oud S arable crop 
E Generally, the model was able to predict seasonal 
2 Ei trends in Rp at most of the sites; however, the model 
& E E occasionally over/underestimated the flux values dur- 
5 EK S 23 3% 8 SG. 2 ing the warm weather in spring and summer. This is 
EE a 2 ist et g particularly evident at the Lincolnshire site, resulting in 
S a high mean difference between modelled and EC- 
F 3 e = a Kee z derived R, (Table 4). Despite using a generic method to 
2e¢g/S 2 & 2ng a estimate Rp from Reco, therefore providing a challenging 
S test for the model, we found no significant mean differ- 
p 2 "S ence between modelled and EC-derived Ru, at three sites 
K 2 BI eo = CH co ch 5 (for a total of four land uses), proving that the model 
> ZO a a an g ; , , 
E 5 adequately simulates soil processes under different 
KR fj land-use systems and climate/soil conditions. 
£ S $ e Low correlation between measurements and model 
a = E S simulations arose predominantly when comparing 
x O g 00 m N oo g model outputs against the IRGA-derived data set; this is 
SS A le ó a oo 5 mainly due to the nature of the measurements (single 
T H data point representing total monthly CO, flux), an 
E £ G aspect not related to the soil processes described in the 
E E a model. However, it is to notice that the IRGA-derived Rp 
& g a a a ar a has been estimated from direct measurements of total 
K aD O ad o a PRES 2 soil respiration and the degree of correlation between 
g E measured and modelled R, is also related to the Rp : Riot 
z cd E ratio adopted. On the other hand, the EC-derived Rp 
E E 3 ^ was estimated from the Reco during daytime, which is a 
By Š E E modelled flux driven by air temperature and other envi- 
E Š 2 z È S 8 © + E ronmental factors. Further model evaluation should be 
£ K based on comparison of the model output with direct 
e 3 Sb measurements of soil Ry fluxes, possibly using automatic 
Z Z £ = Ai E a chambers on soil plots where roots have been excluded. 
zs A S SZ E 5 Kl s This measurement technique would provide continuous 
S S E ; E : S 2 Z 5 R E Ry measurements which would be directly comparable 
= SZ JUZ. sree SEB 3 to the model outputs and therefore would provide a 
E S SE ere 9 Sna Z s| < more accurate evaluation of the performance of the 


© 2015 The Authors. Global Change Biology Bioenergy Published by John Wiley & Sons Ltd., doi: 10.1111/gcbb.12298 


14 M. DONDINI et al. 


model. However, given the very limited input data used 
to run the model and the number of sites/locations used 
for the model evaluation, we conclude that the simula- 
tions are robust and the model adequately simulate soil 
CO, fluxes under five land-use systems. 

Model simulations of N-O and CH; fluxes resulted 
in low correlation and association at most of the study 
sites (Tables 5 and 6), which is expected with such 
low fluxes, and does not represent a failure of the 
model. In fact, the measured N-O and CH, fluxes are 
pooled from sample data points containing outliers 
and extreme variation between sample points in each 
site, which results in a high standard error of the 
measured values. But the N-O and CH: flux simula- 
tions are within the 95% confidence interval of the 
measured values, showing that the model cannot be 
improved to better fit these data and suggesting that 
the lack of correlation between modelled and 
measured values is due to the high variation in the 
measured fluxes, which is a common phenomenon 
verified in many NO (e.g. Oenema et al., 1997; Skiba 
et al., 2013; Cowan et al., 2015) and CH; flux measure- 
ment experiments (Parkin et al., 2012; Savage etal., 
2014). Moreover, if the measured values do not show 
any seasonal trend, a significant correlation with the 
model outputs cannot be obtained (Smith & Smith, 
2007) and low correlation is expected. 

Measured fluxes of CH, were shown to be negligible 
across all land uses and their contribution to the total 
GHG balance, when converted to CO, equivalent, was 
on average <0.2%, except for the Miscanthus field at the 
Aberystwyth site (3% of the total GHG balance). The 
high mean value recorded for Miscanthus in 2012 is 
driven by one replicate with very high CH, production 
and there was large standard error associated with the 
measurements. In general, CH, production or consump- 
tion was negligible also for this field. 

Across all land uses, measured fluxes of NO 
represent a small proportion (<1.5%) of the total GHG 
balance, with the exception of the arable field at the Lin- 
colnshire site and the Miscanthus field at the Aberyst- 
wyth site (6% of the total GHG balance over the 2 years 
measurement period at both fields). Due to technical 
issues and issues regarding access to sites for sampling, 
the data set for the arable and SRC-Willow fields at East 
Grange is missing a substantial number of months, and 
therefore, it was not possible to determine the annual 
GHG balance. 

Despite the very low values of the CH, and NO 
fluxes, and their small contribution to the total GHG 
balance at all experimental sites, both fluxes have been 
modelled adequately on a monthly time-step and no 
improvements can be made to the model with the avail- 
able flux data. 


In this study, all major GHG fluxes from five land-use 
systems were reasonably well estimated using the 
ECOSSE model. The results from this evaluation exer- 
cise show that ECOSSE is robust for simulating GHG 
fluxes from cropland, grassland, SRC-Willow, SRF-Scots 
Pine and Miscanthus (and transitions from the former 
two land uses to the latter three energy crops). This vali- 
dation builds confidence that the model can be used to 
investigate the impacts of land-use transitions spatially 
in the UK and to investigate the effects of converting 
large areas to grow bioenergy crops. 
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